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Abstract 

Results of linear bifurcation and nonlinear analyses 
of the Space Shuttle superlightweight (SLWT) external 
liquid-oxygen (L0 2 ) tank are presented for an important 
end-of-flight loading condition. These results illustrate an 
important type of response mode for thin- walled shells, 
that are subjected to combined mechanical and thermal 
loads, that may be encountered in the design of other liq- 
uid-fuel launch vehicles. Linear bifurcation analyses are 
presented that predict several nearly equal eigenvalues 
that correspond to local buckling modes in the aft dome of 
the L0 2 tank. In contrast, the nonlinear response phenom- 
enon is shown to consist of a short-wavelength bending 
deformation in the aft elliptical dome of the L0 2 tank that 
grows in amplitude in a stable manner with increasing 
load. Imperfection sensitivity analyses are presented that 
show that the presence of several nearly equal eigenvalues 
does not lead to a premature general instability mode for 
the aft dome. For the linear bifurcation and nonlinear anal- 
yses, the results show that accurate predictions of the re- 
sponse of the shell generally require a large-scale, high- 
fidelity finite-element model. Results are also presented 
that show that the SLWT L0 2 tank can support loads in 
excess of approximately 1 .9 times the values of the oper- 
ational loads considered. 

Introduction 

The International Space Station (ISS) is currently 
planned to occupy a 51.6° orbit. Construction of the ISS 
will require the Space Shuttle to deliver a large number of 
payloads to this high-inclination orbit. However, achiev- 
ing this orbit requires that the payload capacity of the or- 
biter be reduced by approximately 10,000 lb. To recover 
part of this lost payload capacity, and to minimize the 
number of Space Shuttle flights needed to build the ISS, 
the National Aeronautics and Space Administration 
(NASA) is developing a new lightweight external fuel 
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tank for the Space Shuttle. This new design, referred to 
as the superlightweight external tank (SLWT), is made 
primarily of an aluminum-lithium alloy and is expected 
to weigh approximately 58,000 lb, which is approximate- 
ly 8,000 lb lighter than the aluminum external tank cur- 
rently in service. This 8,000 lb weight savings translates 
into an 8,000 lb increase in the payload capacity for the 
orbiter. 

An important consideration in the design of the 
SLWT is the nonlinear behavior of its thin-walled re- 
gions that experience compressive or shear stresses, and 
the sensitivity of this behavior to initial geometric imper- 
fections. Small initial geometric imperfections are 
known to cause premature buckling of thin shell struc- 
tures. These effects are very important in the present 
study because local or global buckling of the SLWT shell 
wall could cause the thermal protection system (TPS) to 
separate from the tank, which could cause the vehicle to 
fail. To eliminate shell-wall instabilities from the design, 
accurate predictions of the nonlinear response and im- 
perfection sensitivity of the SLWT are needed. However, 
obtaining this information is a significant task. For exam- 
ple, accurate predictions of the nonlinear response of the 
SLWT have been shown in Ref. 1 to require a large- 
scale, high-fidelity finite-element model to represent the 
complex structural details of the SLWT and a robust non- 
linear shell analysis capability that can predict local and 
general instability buckling modes. 

One component of the SLWT that experiences sig- 
nificant compressive stresses is the liquid-oxygen (L0 2 ) 
tank (see Fig. 1). Prior to launch, the weights of the liq- 
uid-hydrogen (LH 2 ) tank, the L0 2 tank, and the fuel are 
reacted at the solid-rocket-booster attachment points, 
which causes meridional compressive stresses and shear 
stresses that extend into the nose of the SLWT. The non- 
linear behavior of the SLWT L0 2 tank subjected to two 
critical prelaunch loading conditions has been document- 
ed extensively in Ref. 1 . The tank also experiences simi- 
lar compressive and shear stresses during ascent, before 
the two solid rocket boosters (SRBs) are jettisoned. After 
the SRBs are jettisoned, and prior to orbital insertion, the 
L0 2 tank experiences compressive stresses in the aft end 
of the tank instead of the nose region. 

The present paper presents results of linear bifurca- 
tion and nonlinear analyses of the L0 2 tank that were con- 
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ducted at the NASA Langley Research Center. The 
results are for a critical flight loading condition that oc- 
curs 423 seconds into the flight, after the SRBs are jetti- 
soned and prior to orbital insertion (referred to herein as 
the 3g end-of-flight loading condition). First, an overview 
of the SLWT structure and the details of the loading con- 
dition are presented. Then, details of the finite-element 
model, mesh convergence studies, and load simulation 
are summarized. Next, details of the linear bifurcation 
analyses are presented, and then results of nonlinear anal- 
yses for geometrically perfect and imperfect shells are 
presented. The imperfection sensitivity results include a 
discussion of the effects of modal interactions associated 
with several nearly equal eigenvalues. Finally, generic as- 
pects of the finite-element model, analyses, and results 
that may be applicable to the design of future liquid-fuel 
launch vehicles are discussed. 

Overview of the Structure 

The Space Shuttle consists of the orbiter, two SRBs, 
and the external tank (ET), as shown in Fig. 1. The ET 
consists of a L0 2 tank, a LH 2 tank, and an intermediate 
structure called the intertank (Fig. 1). The intertank trans- 
mits the weight of the fuel, the ET structural weight, and 
the orbiter weight to the SRBs prior to launch, and trans- 
mits thrust loads from the SRBs and the orbiter to the ET 
during ascent. The SLWT L0 2 tank is a thin- walled 
monocoque shell that is made primarily of 2195 alumi- 
num-lithium alloy. The L0 2 tank is approximately 49.4 
ft long and has a maximum diameter of approximately 
27.6 ft, as indicated in Fig. 2. The L0 2 tank consists of a 
forward ogive section made from eight gore panels, an aft 
ogive section made from 12 gore panels, a cylindrical bar- 
rel section made from four barrel panels, and an aft ellip- 
tical dome section made from 12 gore panels. The 
coordinate systems used to locate the elements of the LC^ 
tank and the intertank are also shown in Fig. 2. The coor- 
dinates (XT, Y, Z) are typically referred to as the global 
coordinate system of the ET, and axial positions along the 
tank are indicated by the coordinate value of XT in units 
of inches. For example, the location of the junction be- 
tween the forward and aft ogives is indicated by writing 
XT = 536.74 in. Cylindrical coordinates are also used and 
are given by (XT, r, 0), where a positive value of 0 is 
measured from the positive Z-axis toward the positive Y- 
axis, as shown in Fig. 2b. 

The L0 2 tank also has a forward ring frame with a 
"T M -shaped cross section that is referred to herein as the 
T-ring frame, and an aft ring frame with a "Y"-shaped 
cross section that is referred to herein as the Y-ring 
frame. These two ring frames support a baffle assembly 
that prevents the fuel from sloshing during ascent. The 
slosh baffle, a lightweight (approximately 455 lb), thin- 



External tank 

Solid rocket booster (SRB) 


Orbiter 




Forward orbiter 
attachment 
fitting 


Aft orbiter 
attachment 
fitting 


Nose cone 


Forward 
ogive 
section 

Aft ogive 
section 

Cylindrical 
barrel 
section 

Aft elliptical 
dome 
section 

Thrust panel 
SRB Beam 


Aft SRB 
attachment 
point 


Fig, 1 Space Shuttle External Tank components. 

walled structure, is supported by deep, thin- walled rings 
at each end that attach to the forward T-ring and the aft 
Y-ring frames. Other parts of the L0 2 tank include a 
nonstructural nose cone, a forged forward ogive fitting 
and cover plate, an aft spherical dome cap that contains 
the L0 2 suction fitting and a covered manhole, and a vor- 
tex baffle attached to the base of the aft dome cap. The 
L0 2 tank gore and barrel panels are stretch formed, 
chemically milled, and then welded together. The panels 
are fabricated with substantial thickness tailoring to re- 
duce structural weight. The panels are somewhat thicker 
at the welds to form a stiffener-like region that is used as 
a weld land. The primary role of the weld lands is to 
compensate for any reduction in shell-wall strength that 
is caused by welding. Tapering the weld lands in thick- 
ness and width along their length reduces weight and al- 
leviates stress concentrations in the shell that result from 
abrupt changes in thickness. 


The intertank is a right circular cylinder that is 
made from 2090 aluminum-lithium and 7075 aluminum 
alloys and is shown in Fig. 1. The approximately 22.5- 
ft-long intertank has a diameter of approximately 27.5 ft 
and consists of six 45° curved panels that are stiffened 
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Fig. 2 Space Shuttle External L0 2 Tank components 
(values of XT given in inches). 

longitudinally with external hat stiffeners and are re- 
ferred to herein as skin-stringer panels. The intertank 
also has two massive 45° curved panels, referred to as 
thrust panels (see Fig. 1), that are located perpendicular 
to the Y-axis of the intertank and stiffened longitudinally 
with integrally machined external blade stiffeners. 
These eight panels are assembled into the intertank with 
mechanical fasteners and are attached to five large inter- 
nal ring frames, a forward flange, and an aft flange. 
Longitudinal straps (referred to herein as roll ties) sup- 
press lateral-torsional deflection of the ring frames. The 
main central ring frame, two thrust panel longerons, and 
the thrust panels are connected to each end of a tapered 
beam that is referred to herein as the SRB beam (see Fig. 
1). The SRB beam spans the diameter of the intertank 
along the Y-axis and has a maximum depth (in the XT di- 
rection) of approximately 43 in. at its midspan. Forged 
fittings (referred to herein as SRB thrust fittings) that are 
incapable of transmitting moments are fastened to the 
ends of the SRB beam. The primary role of the thrust 
panels is to diffuse the large axial loads introduced by the 
SRBs into the intertank and then into the L0 2 tank shell 


wall. The SRB beam compensates for the eccentricity of 
the concentrated loads introduced by the SRBs. The 
SRB beam also supports loads that are normal to the in- 
tertank (parallel to the SRB beam) at the SRB attachment 
points. The intertank also has a 46-in.-high by 52-in.- 
wide frame-reinforced nonstructural access door located 
along the cylinder generator at approximately 6 = 146°. 

Critical Loading Condition 

The 3g end-of-flight loading condition was identi- 
fied by the members of the SWLT team at the NASA 
Marshall Space Flight Center and the Lockheed Martin 
Manned Space Systems Company as a critical loading 
condition that could cause buckling of the aft dome of the 
L0 2 tank. This critical loading condition corresponds to 
a partially filled L0 2 tank aft dome and an acceleration 
field given by a x = -3.040g, ay = 0.005g, and a z = 
-0.725g, where g is the magnitude of the gravitational ac- 
celeration. The X subscript in the first acceleration cor- 
responds to the XT direction and is used in the present 
paper for convenience. Details of this loading condition 
are shown in Figs. 3 and 4. The loads consist of the iner- 
tial loads of the structural mass and the L0 2 mass, the ul- 
lage pressure exerted on the shell wall (see Fig. 4), the 
LH 2 tank interface force and moment, and the thermal 
load associated with the cryogenic fuel and aerodynamic 
heating. No SRB forces or aerodynamic pressure loads 
are present for this loading condition because the SRBs 
have been jettisoned by this point in the flight profile and 
the orbiter and external tank are essentially out of the 



Fig. 3 Loads at 423 seconds into flight. 
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dense part of the Earth’s atmosphere. The inertial loads 
of the structural mass (magnitudes) are given by m*. a x = 
72.38 kips, m s a Y = 0.12 kips, and m s a z = 17.26 kips. 
Similarly, the inertial loads of the L0 2 mass are given by 
m L02 a x = 588.67 kips, m L02 a Y - 0.99 kips, and m L02 
a z = 143.18 kips (see Fig. 3). The interface force and 
moment between the intertank and the LH 2 tank are giv- 
en by F = - 661.0501+ 1.1 lOj -160.440k kips and M = 
-6.040i- 5 1, 428 .720j -329.220k in-kips, where i,j,and 
k are standard orthonormal base vectors associated with 
the XT, Y, and Z axes, respectively. 

As the ET is accelerated, the L0 2 mass exerts pres- 
sure on the L0 2 tank shell wall. The resultant force of 
this pressure distribution corresponds to the inertial loads 
of the L0 2 that are given in Fig. 3. This pressure distri- 
bution is approximated by the situation depicted in Fig. 
4 in which the L0 2 and its container are subjected to rig- 
id-body motion and flow effects are neglected. This ap- 
proximate pressure distribution is given by 


p(x,0) = p u for x s x f ± r(x)|j“ sin0 + ~ cosej and 

p(x,0) = p u ±YixJg X ( x± *f) + tOO^jr sin6 + cose 


where x is a local 


for x>x f ±i{xM^sin0 + ^cos0 
1 \ a x a x f 

axial coordinate that is measured from XT = 371 .00 in., 

x f = 498 in. is the local x-coordinate of the fill level (XT 

= 869.00 in.), and r(x) is the horizontal or polar radius of 

the tank. The ullage pressure is given by p u = 17.0 psig 

and the specific weight of the L0 2 that was used in the 

present study is given by Ylo2 = 0.04123 lb/in . The 


local z -axis shown in Fig. 4 corresponds to the axis of 
maximum inclination of the free surface of the L0 2 and 

is given by 0 = 180° + ip, where = Tan ±1 (a Y /a z ) = 

- 0.40°. Similarly, the free surface of the L0 2 has an 




Fig. 4 L0 2 pressure distribution at 423 seconds 
into flight. 


inclination angle in the x-z plane that is given by 

5 = Tan*|s/ a Y + a| / | a x|j = 13.41°. 

The temperature distribution that was used in the 
present study as an approximation associated with the 
cryogenic fuel and aerodynamic heating is shown in Fig. 
5. The temperature distribution shown in this figure is 
axisymmetric and varies along the surface meridians in a 
piecewise-linear manner. The hottest place on the L0 2 
tank and the intertank is the tip of the ogive (275 °F), and 
the coldest is at the bottom region of the aft dome (-297 
°F). The nominal ambient temperature of the L0 2 tank 
and the intertank prior to fueling is 50 °F. This temper- 
ature was used in the present study as the temperature at 
which thermal stresses in the L0 2 tank are absent. 

Analysis Code and Finite-Element Modeling 

The results of the linear bifurcation buckling and 
nonlinear analyses were obtained with the STructural 
Analysis of General Shells (STAGS) nonlinear structural 
analysis code for general shells. 2 The finite-element mod- 
els of the SLWT tank that were used in the present study 
are very complex and include many structural details and 
the skin thickness variations or tailoring used to reduce 
structural weight. A detailed description of these models 
is presented in Ref. 1 . STAGS was chosen for analyzing 
the SLWT tank because of its robust state-of-the-art non- 
linear-equation solution algorithms and its general user- 
input capability that is convenient for modeling branched 
shells typically used for launch vehicles. In particular, 
STAGS uses both the full and modified Newton methods 
to obtain an accurate nonlinear solution, and large rota- 
tions in the shell are represented by a co-rotational algo- 
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Fig. 5 Axisymmetric temperature profile at 423 
seconds into flight (values of XT are given in in.). 


rithm at the element level. The Riks arc-length projection 
method is used to continue a solution past limit points. 
STAGS permits complex geometries, loading conditions, 
and initial geometric imperfections to be modeled in a di- 
rect manner by the use of user-written subroutines that are 
essentially independent of the mesh discretization. For 
example, these user- written subroutines allow the user to 
define reference surface geometries; tapered shell walls 
and stiffener cross-sections; and complex nodal force, 
temperature, and pressure distribution functions in a di- 
rect manner, using the FORTRAN computer language. 
This feature greatly simplified the definition of the finite- 
element models and the mesh convergence studies con- 
ducted in the present study. A description of how the fea- 
tures of STAGS were used in the present study to model 
the SLWT L0 2 tank and intertank, and details of how ap- 
plied prelaunch loads were simulated are presented in 
Ref. 1 . 

The basic approach used in the present study to 
simulate the actual 3g end-of-flight loading condition is 
to apply all accelerations and associated inertial and L0 2 
pressure loads illustrated in Figs. 3 and 4 to the model. 
The thermal load was applied by introducing the axisym- 
metric temperature distribution shown in Fig. 5 as a tem- 
perature change from a nominal initial uniform 
temperature of 50°F. The bottom of the intertank, where 
the LH 2 tank interface force and moment act, was re- 
strained so that the LH 2 tank interface force and moment 
become reactions, and rigid body motion is eliminated. 
Next, the applied loads were separated into two groups. 
The first group contains the pressure that acts on the shell 
wall because of the acceleration of the L0 2 mass, the 
structural-mass inertial loads, and the inertial line loads 
that represent the acceleration of the slosh baffle mass 


that is located inside of the barrel section of the L0 2 tank. 
This group of loads is treated as the primary source of de- 
stabilizing compressive stresses in the L0 2 tank that may 
occur at load levels greater than the corresponding oper- 
ational load level. The second group of loads consists of 
the thermal load and the L0 2 tank ullage pressure. The 
loads in the second group are constant in value, are part 
of the operational loads, and are considered to be passive 
loads when determining the stability margin of safety of 
the L0 2 tank. In performing linear bifurcation buckling 
and nonlinear analyses with STAGS , two load factors, p a 
and pfc, were assigned to the first (active) and second 
(passive) load groups, respectively. Values of p a = pb 
= 1 correspond to the loading condition that was men- 
tioned previously as the operational load level. 

In modeling the SLWT L0 2 tank and the intertank, 
several assumptions were made to simplify the finite-el- 
ement models. Limited parametric studies were con- 
ducted to determine the adequacy of the assumptions and 
simplifications. For each case in these studies, the mod- 
eling assumptions used to simplify the finite-element 
models were found to be acceptable for analyzing the 
nonlinear behavior of the SLWT L0 2 tank. In addition, 
the finite-element modeling approach that was used for 
the SLWT L0 2 tank was also applied for the analysis of 
two full-scale structural tests that were conducted at the 
NASA George C. Marshall Spaceflight Center on the 
original standard weight ET during the development pro- 
gram of the original Space Shuttle ET. The analytical re- 
sults for these two test articles, which buckled 
unexpectedly during the tests, indicate that the finite-el- 
ement modeling approach that is used in the present 
study is adequate for representing the nonlinear behavior 
of the SLWT L0 2 tank (Ref. 1 ). 

Results and Discussion 

Three different finite element models were used in 
the present study for analysis of the L0 2 tank subjected 
to the 3g end-of-flight loading condition. As a first step 
toward identifying an adequate model with as few de- 
grees of freedom as necessary, linear bifurcation buck- 
ling analyses were conducted. The passive loads 
associated with load factor Pb were applied to the 
STAGS models as a linear prebuckling stress state (pb = 
1) and the active (destabilizing) loads associated with 
load factor p a were used to obtain the minimum eigen- 
value. The models that were investigated had 76,900, 
103,400, and 158,000 degrees of freedom. The model 
that was identified as adequate for predicting the linear 
bifurcation buckling behavior is shown in Fig. 6 and cor- 
responds to 158,000 degrees of freedom. 

The first linear bifurcation mode (referred to herein 
as the linear bifurcation buckling mode) for the geometri- 
cally perfect shell is shown in Fig. 7 for the STAGS mod- 
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(a) Mesh with intertank shown. 



Fig. 6 Finite element mesh (158,000 degrees of freedom). 

el with 158,000 degrees of freedom. The eigenvalue for 
this model corresponds to active loads that are approxi- 
mately 1 .526 times the magnitude of the corresponding 
operational loads shown in Figs. 3 and 4. This buckling 
mode is a localized, short-wavelength wrinkle in the for- 
ward part of the aft dome that extends from approximate- 
ly 0 = 98° to 132° with its large amplitude region centered 
approximately on 0 = 1 12.5°. The second through fourth 
linear bifurcation modes are also localized, short- wave- 
length modes, similar to the mode shown in Fig. 7, with 
eigenvalues equal to 1 .528, 1 .540, and 1 ,541 , respective- 
ly, which are all less than 1% higher than the lowest 
eigenvalue. The second linear bifurcation mode, which is 
nearly identical to the first mode, is also in the forward 
part of the aft dome with its large amplitude region cen- 
tered very close to 0 = 112.5°. The third and fourth 
modes are also nearly identical and located in the forward 
part of the aft dome, but their large amplitude regions are 
centered approximately on 0 = 247.5°. 

The nearly equal values of the first four eigenvalues, 
the short wavelength of the linear bifurcation modes, and 
the large difference in location between the first two and 
last two linear bifurcation modes led to the dense mesh re- 
finement of the aft dome that is shown in Figs. 6 and 7. A 
uniform circumferential mesh refinement of the dome 
was used in order not to bias the location of the response, 
and was facilitated by the use of the five-node and seven- 
node rectangular transition elements available in 
STAGS. 3 The 103,400- and 158,000-degree-of-freedom 
models have the same general mesh arrangement shown 
in Fig. 6, but the level of refinement of the forward part of 
the aft dome shown in Figs. 6 and 7 for the 158,000-de- 



Fig. 7 Linear biforcation buckling mode (158,000 degrees 
of freedom; p a = 1526 and Pt, = 1.0). 

gree-of-freedom model is essentially twice that of the 
1 03 ,400-degree-of-freedom model. The lowest eigenval- 
ues for the 103,400- and 15 8, 000-degree -of- freedom 
models are given by p a = 1 .536 and p a = 1 .526, respec- 
tively. The smoothness of the buckling mode shown in 
Fig. 7 and the one-percent difference in the eigenvalues, 
suggest that the 158,000-degree-of-freedom model is ad- 
equate for representing the linear bifurcation behavior of 
the LC >2 tank for this loading condition. The 103,400- and 
158,000-degree-of-freedom models were also used to ob- 
tain nonlinear solutions for geometrically perfect and im- 
perfect shells. These solutions, which are in good 
agreement, indicate that the 158,000-degree of freedom 
model adequately represents the nonlinear behavior of the 
L0 2 tank for this loading condition. Thus, all subsequent 
results presented in this section were obtained with the 
1 5 8, 000-degree -of-freedom model . 

The meridional stress resultant distribution (given 
in units of Ib/in.) in the aft dome, that was obtained from 
nonlinear analyses, is tensile and is shown in Fig. 8 for 
values of p a = p*, = 1 and for p a = 1 .5 and p b = 1. The 
lighter meridional bands of shading shown in the figure 
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Fig. 8 Meridional stress resultants in aft dome (lb/in.). 

correspond to lower values of the meridional stress re- 
sultants along the weld lands. The meridional tension 
shown in the figure primarily resists the axial accelera- 
tion of the L0 2 mass and is responsible for the shortness 
of the wavelength of the linear bifurcation buckling 
mode. The circumferential stress resultants in the aft 
dome that generate the linear bifurcation modes are 
shown in Fig. 9 for values of p a = Pb = 1 and for p a = 1 .5 
and Pb = 1 ■ These results show a region of circumferen- 
tial compression (90° 270°) that is also in the for- 

ward part of the aft dome, with the highest values of 
compression approximately near 0 = 180°. The average 
shell wall thickness in and around this region is 0.090 in. 

An important concern that arose during the course 
of the present study is the possiblity of high sensitivity to 
initial geometric imperfections, that may be caused by the 


(b) p a = 1.5 and p b = 1.0. 

Fig. 9 Circumferential stress resultants in aft dome 

(lb/in.). 

presence of four closely spaced eigenvalues of the linear 
bifurcation modes. This sensitivity could lead to a prema- 
ture collapse mode of the aft dome. This concern led to 
the use of an imperfection shape that is in the form of a 
linear combination of the first four linear bifurcation 
modes, described previously in the present paper, in the 
nonlinear analyses of the imperfect shell. With this four- 
mode imperfection each mode possesses an equal degree 
of participation. Mathematically, this set of modes can be 
viewed loosely as a basis for an "isotropic imperfection 
space,” similar to a basis of a vector space. The first four 
modes were selected because they represent deformation 
states that the structure has an intrinsic affinity to deform 
into, provided that there are no substantial nonlinear pre- 
buckling effects present. That is, in the absence of sub- 




stantial nonlinear prebuckling effects, the linear 
bifurcation eigenvalues represent when (at what load lev- 
el) and where (what configurations) strong interactions 
between compressive membrane stresses and normal dis- 
placements are likely to be present. In addition, the four- 
mode imperfection was selected because of the statement 
given by Bushnell 4 that suggests that premature failure of 
shell structures that exhibit a short-wavelength response 
can be activated or "triggered” by imperfections with a 
similar short-wavelength shape. Thus, the four- mode im- 
perfection is expected to represent adequately a preferred 
direction of departure from the primary equilibrium path 
should the structure have a tendency to do so. 

Results are presented in Figs. 10 and 1 1 that show 
the nonlinear deformations that were obtained from 
STAGS analyses of a geometrically perfect shell and a 
geometrically imperfect shell with an imperfection-am- 
plitude-to- wall-thickness ratio AJt = 1 , respectively. The 
thickness t in the ratio AJt is the average wall thickness 
of the top part of the aft dome (t = 0.090 in.). The results 
shown in Fig. 1 1 were obtained by using the four-mode 
imperfection that was described previously in the present 
paper. For this imperfection, the imperfection shape is 
given by -0.707 times the summation of the first four 
eigenvectors. The negative of the linear combination was 
used as the imperfection shape because it was found to 
provide a stronger nonlinear interaction with the com- 
pressive stresses in the shell wall than the positive-valued 
linear combination. The results shown in these two fig- 



Fig. 10 Nondimensional normal displacement w/t of 
geometrically perfect aft dome (t = 0.090 in.; XT = 885 
in.). 
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Fig. 11 Nondimensional normal displacement w/t of 
geometrically imperfect aft dome. (Imperfection- 
amplitude-to- wall-thickness ratio AJt = 1; t = 0.090 in.; 
XT = 885 in.). 

ures are for nonlinear solutions that were obtained by in- 
creasing the load factors p a and p^ simultaneously to a 
value of one, and then holding p*, constant while increas- 
ing the magnitude of the load factor p a . The nondimen- 
sional normal displacements w/t along the circumference 
of the aft dome shell wall at XT = 885 in. and from 0 = 
80° to 1 80° are represented by the solid lines in the fig- 
ures for values of the load factor p a approximately equal 
to 1.0, 1.5, and 1.8. Moreover, the line given by 0 = 180° 
in each figure is essentially a line of reflective symmetry 
for the nonlinear-response curves. Overall, negative val- 
ues of the normal displacements are indicated by the left- 
hand ordinate for these three load factors. These results 
are negative because of the contraction of the aft dome 
that is caused primarily by the LO 2 thermal load (shrink- 
age). The linear bifurcation buckling mode and half of 
the four-mode imperfection are represented by the dashed 
lines in Figs. 10 and 11, respectively, with normalized 
amplitudes given by the right-hand ordinate of the fig- 
ures. The third linear bifurcation mode is obtained by re- 
flection of the dashed line in Fig. 10 about the line 0 = 
180°. Similarly, the second half of the four-mode imper- 
fection is obtained by reflection of the dashed line in Fig. 
1 1 about the line 0 = 180°. The linear bifurcation buck- 
ling mode and the half of the four-mode imperfection 
shown in the figures are very similar and are included in 
the figures to give an indication of how the imperfection 
shape influences the nonlinear solution. 

The solid lines shown in the Figs. 10 and 11 indi- 
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cate a short-wavelength bending response in the circum- 
ferentially compressed part of the aft dome that extends 
from approximately 0 =95° to 180°. Although they are 
not shown in the figures, the deformations are symmetric 
about the right-hand axis and actually extend to approx- 
imately 0 = 265°. The pattern of the nonlinear deforma- 
tions are substantially different from the linear 
bifurcation buckling mode and the four-mode imperfec- 
tion, and encompass a much larger portion of the dome 
circumference. The overall increasing slope trend of the 
solid lines in the figures (seen by neglecting the undula- 
tions in the curves) is due to the inward displacement of 
the shell wall that is caused by the increase in pressure 
and temperature change as p a and Pb are simultaneous- 
ly increased to a value of one in the nonlinear analysis. 
This effect is not accounted for in a linear bifurcation 
buckling analysis and, as a result, the dashed lines do not 
exhibit a similar overall increasing slope trend. 

The results presented in Figs. 10 and 1 1 predict a 
stable, but highly wrinkled, nonlinear response at a load 
level greater than the load level predicted by the linear 
bifurcation buckling analysis (p a = 1 .526). At the opera- 
tional load level given by p a = 1 , the results predict sig- 
nificant nonlinear prebuckling deformations at values 
close to 0 = 105°, 135°, and 165° which are the locations 
of the meridional weld lands in that part of the aft dome. 
As the load is increased, substantial bending deforma- 
tions (indicated by the waviness of the curves) develop 
and grow in the shell wall between the weld lands, which 
reduces the apparent circumferential stiffness of the aft 
dome. The largest bending deformations shown in the 
figures are predicted to occur at the locations of the me- 
ridional weld lands and are thought to be slightly over- 
predicted because the weld lands are represented by 
discrete stringers (concentrated stiffnesses). Although it 
is difficult to see in the figures, these results also indicate 
that a geometrical imperfection in the shape of the first 
four linear bifurcation modes and with a small negative 
amplitude does amplify the severity of the bending de- 
formations (e.g., compare the magnitudes of Aw/t shown 
in the figures) and causes the growth of the bending de- 
formations to start at slightly lower load levels. 

The reduction in the apparent circumferential stiff- 
ness of the aft dome is shown more explicitly in Figs. 12 
and 13 for the geometrically perfect and imperfect shells, 
respectively. In these figures, the intensities of the bend- 
ing deformations along XT = 885 in., that are shown in 
Figs. 10 and 1 1 , are given as a function of the load factor 
p a for several circumferential locations. The intensity of 
the bending deformations is given by the nondimension- 
al quantity Aw/t, which is the normalized amplitude of 
an undulation in the curves at a specific value of 0. For 
example, Aw/t is shown in Figs. 10 and 11 for approxi- 
mately 0 = 110°. Four curves are presented in Fig. 12 


6 = 180 ° 121 ° 



Fig. 12 Local nondimensional normal displacement 
amplitudes Aw/t of geometrically perfect aft dome at 
XT = 885 in. (t = 0.090 in.). 


for the geometrically perfect shell that correspond to val- 
ues of 0 = 110°(filled circles), 121°(unfilled triangles), 
169°(unfilled squares), and 180°(filled triangles). The 
curves for 0 =110° and 169° correspond to the severe 
bending deformations next to two of the weld lands, and 
the other two curves correspond to locations that are ap- 
proximately midway between weld lands. Similarly, the 
curves for 0 = 1 10° and 180° that are shown in Fig. 12 
are also shown in Fig. 13 (solid curves) along with the 
corresponding curves for the imperfect shell with A/t = 
1 (dashed curves). The horizontal dashed lines shown in 
Figs. 12 and 13 represent the linear bifurcation buckling 
load level (p a = 1.526). 

The results shown in Figs. 12 and 13 indicate that 



Fig. 13 Local nondimensional normal displacement 
amplitudes Aw/t of geometrically perfect and imperfect 
aft domes at XT = 885 in. (t = 0.090 in.). 
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the bending deformations near the meridional weld lands 
dominate the response of the aft dome; that is, Aw/t 
grows at the highest rate at 6 = 110° and 169°. More- 
over, the results predict that this effect is much more pro- 
nounced for the imperfect shell with A/t = 1 . The curves 
for 0 = 1 10° were included in Fig. 13 because they show 
the greatest effect of the four-mode imperfection. In 
contrast, the effect of the imperfection is shown to be be- 
nign for the point between weld lands at 0 = 180°. Cor- 
responding results were obtained for 0 =121° and 169°. 
These results show that the effect of the imperfection at 
0 =121° (between weld lands) was slightly larger than 
that shown in Fig. 1 3 for 0 =1 80°, but much smaller than 
that shown for 0 =110° (next to a weld land). Also, the 
results show that the effect of the imperfection at 0 = 
169° (next to a weld land) is benign. 

All the results shown in Figs. 12 and 13 show a 
monotonically increasing nonlinear response and predict 
that the shell can support loads greater than the critical 
buckling load predicted by a linear bifurcation buckling 
analysis. As Aw/t increases, the apparent circumferen- 
tial stiffness decreases, and as a result, the positive-val- 
ued constant of proportionality between an increment in 
load and the corresponding increment in displacement 
amplitude decreases. This trend is manifested by the re- 
duction in slope of the load versus displacement ampli- 
tude curves. This type of response is similar to the 
response presented for the prelaunch loading condition 
with full L0 2 and LH 2 tanks, 1 and to the response report- 
ed by Stevens, Starnes, and Almroth 5 for cylindrical 
shells subjected to combined internal pressure and a pure 
bending moment. The results in Ref. 5 indicate that the 
amplitude of the short-wavelength deflection grows rap- 
idly as the load increases and approaches a critical value. 
At the critical value of the load, the load-deflection re- 
sponse curve approaches a horizontal tangent that corre- 
sponds to a local collapse mode of the cylinder. 
Mathematically, the horizontal tangent indicates that un- 
bounded growth of the displacement occurs for an infin- 
itesimal increase in the load. It is expected that the 
curves shown in Figs. 12 and 13 would approach a hori- 
zontal tangent as Aw/t increases until a redistribution in 
load occurs within the aft dome. As a horizontal tangent 
in a load versus displacement amplitude curve is ap- 
proached, the region of the shell containing the bending 
deformations becomes incapable of supporting addition- 
al load, and the compressive load is redistributed to an- 
other portion of the aft dome. If other parts of the dome 
cannot support the redistributed compressive load or if 
excessive yielding occurs, the shell will collapse. 

The results presented in Figs. 10 through 13 indi- 
cate that large local bending deformations may occur in 
the shell wall for loads that are much smaller than the lo- 
cal collapse load and may cause the thermal protection 


system (TPS) to debond from the shell wall and fail. 
This mode of failure is of great importance in the design 
of the TPS for contemporary space vehicles. The results 
presented in Fig. 14 give approximate estimates of the lo- 
cal radius of curvature p for the bending deformation in 
the aft dome at XT = 885 in. and 0 = 109° (next to a weld 
land). The local radius of curvature p was calculated 
with STAGS as close as possible to the crests of the de- 
formations patterns where the slope w ; = 0. At this loca- 
tion, the radius of curvature is given by p = |w /f | ±] , 
where w is the normal displacement and the prime marks 
in the equation denote differentiation with respect to a lo- 
cal Cartesian surface coordinate. The filled circles shown 
in the figure correspond to results for the geometrically 
perfect shell, and the unfilled circles correspond to re- 
sults for the geometrically imperfect shell with imperfec- 
tion-amplitude-to-wall-thickness ratios of A/t = 1 . The 
results in Fig. 14 demonstrate that the geometric imper- 
fection amplitude has a significant influence on the local 
radius of curvature of the shell wall. For example, if a 
given thermal protection system is known to debond 
from the shell wall at a value of p = 100 in., the maxi- 
mum load factor is reduced from a value of approximate- 
ly 1 .6 for the geometrically perfect shell to 1 .3 for the 
geometrically imperfect shell with A/t = 1 . However, a 
large amplitude, short- wavelength imperfection of this 
type is very unlikely to be present in a piece of high-pre- 
cision flight hardware. Thus, the radius-of-curvature re- 
sults for the perfect shell are much more practical than 
the corresponding results for the imperfect shell since 
such large imperfections would most likely be identified 
during inspection of the shell. 

Despite the concerns about acute imperfection sen- 
sitivity that is sometimes associated with the presence of 
several nearly equal linear bifurcation eigenvalues, the re- 



Fig. 14 Local radii of curvature p of bending 
deformation at XT = 885 in. and 6 = 109° for 
geometrically perfect and imperfect aft domes 
(t = 0.090 in.). 
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suits presented in Figs. 10 through 13 indicate a stable 
nonlinear response for the imperfect shell, which has a 
relatively large imperfection amplitude. Moreover, the 
results indicate the presence of substantial nonlinear pre- 
buckling deformations that appear to diminish the signif- 
icance of the closely spaced linear bifurcation 
eigenvalues. This trend is consistent with the qualitative 
guidelines for imperfection sensitivity given in Ref. 4. 
The trend can be explained by noting that, as the nonlin- 
ear prebuckling deformations grow with increasing load, 
they represent a perturbation of the shell geometry and a 
local change in the stress state that is not represented 
physically by the linear bifurcation eigenvalues. The sta- 
bility of the wrinkle-like deformation state is at least par- 
tially attributed to the fact that the compression region is 
a local region that is more likely to cause a benign internal 
load redistribution in the presence of nonlinear prebuck- 
ling deformations than a sudden mode change or collapse. 
This explanation is consistent with qualitative guidelines 
for imperfection sensitivity that are presented in Ref. 4 
and was confirmed in the analyses by performing bifurca- 
tion analyses from the nonlinear configurations at several 
of the higher load levels. The behavior of the aft dome of 
the LO 2 tank is significantly different from that of a com- 
pression-loaded cylinder or an externally pressurized 
sphere which exhibit several nearly equal, or a multiplic- 
ity of, linear bifurcation eigenvalues. A major difference 
is that the region of compression of the aft dome does not 
fully envelop the shell, unlike the compression-loaded 
cylinder or the externally pressurized sphere. This differ- 
ence facilitates load redistribution in the aft dome without 
shell collapse or a mode change . 

Concluding Remarks 

Linear bifurcation and nonlinear analyses of the 
Space Shuttle superlightweight (SLWT) liquid-oxygen 
(LO 2 ) tank have been presented. The loading details for 
an important end-of-flight-loading condition have been 
described and the analytical method used to simulate the 
loading condition has been discussed. Results have been 
presented herein that were obtained from complex, large- 
scale finite-element models of a portion of the Space 
Shuttle SLWT. These results illustrate an important type 
of nonlinear response for thin- walled shells, that are sub- 
jected to combined mechanical and thermal loads, that 
may be encountered in the design of other liquid-fuel 
launch vehicles. In addition, the results indicate that 
large-scale, high-fidelity finite-element models are gen- 
erally required to predict accurately the linear bifurcation 
and nonlinear responses. 

For the end-of-flight-loading condition, linear bi- 
furcation analyses yielded several nearly equal eigenval- 
ues that correspond to local buckling modes. However, 


the nonlinear analyses yielded a response that is charac- 
terized by a short-wavelength bending deformation that 
grows in amplitude in a stable manner with increasing 
load. These local bending deformations appear in the aft 
elliptical dome of the L0 2 tank and do not lead to a gen- 
eral instability mode, but may cause failure of the ther- 
mal protection system (TPS). Because of this concern, 
results have been presented that can be used to determine 
the load level at which a TPS failure is likely to occur. 
Imperfection sensitivity analyses have also been present- 
ed that indicate that the aft dome does not exhibit a non- 
linear collapse mode associated with the interaction of 
nearly equal linear bifurcation modes, for load levels be- 
low approximately 1 .9 times the operational load level . 
However, the results do predict that the severity of some 
of the local deformations is significantly affected by the 
localized initial geometric imperfection. 
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